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1. INTRODUCTION 


Depending on the nature of certain problems, geometers found it expedient to 
look beyond the strict Cartesian systems of coordinates. 

Lamp's investigations mark a brilliant achievement in the history of curvi- 
linear systems, hi addition to his contributicms to differentfal geometry, his now 
classical reference [Lam^, 1837] introduced for the first time the idea of curvilinear 
elliptic coordinates (today the name ''ellipsoidal" is generally accepted). The same 
basic theory has been used by geodesists up to recent times [Molodenskii et al. , I960]. 

After Lame, scores of mathematicians became interested in the subject, 
improving and generalizing the methods. Among them, (Darboux, 1898] some sixty 
years later, wrote several papers originating the concept of moving frames, which 
later was universalized by [Cartan, 1935]. 

In modern times, mathematics has progressed toward the maximum degree 
of generalization. With the advent of Absolute Differential Geometry, the study of 
such abstract topics as m -dimensional manifolds in n-space and the concepts of ten- 
sor and differential forms have drastically revolutionized the field. 

Even in geodesy some pioneering steps intheMe areas have been taken [Marussi, 
1949], [Hotine, 1969], [Grafarend, 1975], culminating in the periodic celebration of 
the Hotine Symposiums on Mathematical Geodesy. 

The present report benefits from some of the above methods in order to develop 
differential transformations between Cartesian and curvilinear orthogonal coordinates. 
However, only matrix algebra is used for the presentation of the basic concepts. The 
fact that second order Cartesian tensors reduce to 3 x 3 matrices frequently is over- 
looked. 

After defining in Chapter 2 the reference systems used in this work. Chapter 3 
introduces the rotation (R), "metric" (H) and Jacobian (J) matrices of the transforma- 
tions between Cartesian and curvilinear coordinate systems. A value of R as a function 
of H and J is presented. Likewise an anal 3 rtical expression for J^as a function of H^ 
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and R is obtained. Subsequently, in Chapter 4, emphasis is placed on showing that 
the diRerential equations published in the English translation of [Molodenskii et al. , 
1960] are equivalent to conventional similarity transformations. This dissipates the 
confusion created recently by some authors [Badekas, 1969], {Krakiwsky and Thompson, 
1974] who credited [Molodenskii et al. , 1962] with a model they never wrote. A 
discussion of scaling methods follows. 

Chapter 5 introduces ellipsoidal coordinates, to which the general theory 
developed in Chapter 3 is applied. Finally, differential transformations between 
ellipsoidal and geodetic coordinates arc established. 
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2. REFERENCE COORDINATE SYSTEMS 


The principal problem of geodesy may be stated as follows (Hirvonen, I960]: 
"Find the space coordinates of any point P at the physical surface S of the earth 
when a sufficient number of geodetic operations have been carried out along S. " 

Therefore, in order to know the position of P, the definition of an appro- 
priate system of coordinates is of primary importance. 

Due to the nature of the rotational motions of the earth and to other geodynamic 
problems, a rigorously defined system of the accuracy of our current observational 
capabilities, is not presently available. A recent colloquium organized by the lAU 
(International Astronomical Union) in Toru^, Poland, was the first attempt to coor- 
dinate the work of different groups in the international scientific community for the 
future definition and selection of reliable reference frames (Kotaczek and Weiffen- 
baeh, 1974]. 

In the present report, only those earth fixed coordinate ^sterns (Terrestrial 
Systems) which are commonly used in geodesy will be described. The reader is 
assumed to be familiar with other celestial systems used frequently in astronomy 
and conveniently defined, for example, in [Mueller, 1969]. With regard to the 
c^namically defined coordinate systems, generally best suited for geophysical prob- 
lems, see the description in [Munk and MacDonalc^ I960]. 

In the first place, a broad division between Cartesian and curvilinear systems 
may be made. Due to the nature of the basic reference surface in geodetic problems, 
sometimes it is convenient <.>se curvilinear coordinates instead of spatial rectangular 
coordinates. This is especially true when the ellipsoid is used as the basic reference. 

hi the following sections the coordinate systems used in this report and their 
notation will be presented in order to avoid any possible confusion. 



2.1 


Quasi-Geocentric Cartesian Systems 


(x, y, z) H "Geographic ** or Mean Terrestrial System 

Origin: Close to the geocenter (center of mass of the earth, including 

the atmosphere) 

z axis: Directed toward the CIO (Conventional International Origin) as 

defined by the IPMS (International Polar Motion Service) and 
the BIH (Bureau International de I'Heure). 

X axis: Passes through the point of zero longitude as defined by the 1968 

BUI system {Guinot et al. , 1971] . 

y axis: Forms a right-handed coordinate system with the x and z axes. 

A redefinition of this system is plausible in the future [see Kdaczck and 
Weiffenbach, 1974, pp, 34-37]. 

World Systems . These are systems defined by particular satellite solu- 
tions accomplished by different organizations. Two wide categories can 
be mentioned: 

Dynamic Solutions (Geocentric) 

(X, y, z)g3f(. Goddard Space Fli^t Center 
(x, y, z) 3 *f, - Smithsonian Astrophysical Observatory 

(x, y, z)n„l = Naval Weapons Suriace Center 
Geometric Solutions (Non-Geocentric) 

(x, y, z)nos National Ccean Survey 
(x, y, z)o 3 u = Ghio State University 

For a complete description of the different published solutions and their 
corresponding references, consult [Mueller, 1975]. 

Some as yet unexplained differences between the orientation of the world 
systems with respect to the geographic system are reported In [Mueller 
et al. , 197.3]. 
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Fig. 2.1 "Geographic," Geodetic and Local Cartesian Systems 



Fig. 2.2 Orthonormal Bases (I, J, k) and (oi, ea» e.-)) 


(u, V, w) = Geodetic Systems (one for each particular local datum) 

Origin: The center of the reference ellipsoid used for defining the 

datum in question. 

w axis: Coincides with the semiminor axis b of the reference ellipsoid, 

u axis: Passes through the point (X = 0, <p-0) 

V axis: Forms a right-handed frame with u and w axes. 

Errors in the deflections of the vertical adopted at the datum origin, in 
the observed astronomic latitude and longitude, and the adoption of 
improper parameters of the referenced ellipsoid shift the origin of this 
system from the geocenter by amounts 6u, 6v, 6w. 

The improper application of the Laplace condition and errors in the 
astronomic a: imuth introduce non-parallelism between the geodetic 
and geographic systems. The relationship is established through the 
rotations 6 c, 6o). See Fig. 2.1. Examples of this type of system 
defined through the datum coordinates are: 

(u, V, w)n/^o - North American Datum 
(u, V, w)guo ^ European Datum 




( X , i8 , u ) = Curviliaear Ellipsoidal (Rotational) Coordinates (as defined in 
[Heiskanen and Moritz, 1967] ) 

X: The same as above (i.e. geodetic longitude ellipsoidal 

longitude). 

/3 : Ellipsoidal or reduced latitude (see Fig. 5.1). 

u : Semiminor axis of the confocal ellipsoid throu^ P. 

(A, H) = Natural or Astronomical Coordinates 

This curvilinear coordinate system refers to the instantaneous terres- 
trial system. In this report only the reduced astronomic coordinates 
will be used. Consult [Mueller, 1969] for the corresponding definitions. 

: Reduced astronomic longitude 

0* : Reduced astronomic latitude 


Frames of Reference 


(tj» C» O - Local Geodetic Frame 

Origin: The point P(X, (p, h). In the case when P is on the 

earth surface, the local coordinate system will be called 
topocentric . 

C axis: Normal through P to the reference ellipsoid. The positive 

sign in the outward direction. 

r) axis: Normal to Z ^^<1 the geodetic meridian plane (vfaen h = 0, 

tangent to the geodetic parallel of P. ) Positive in the direc- 
tion of increasing X. 

Z axis: Perpendicular to t? and Z forming a right-hande- system 

(when h = 0, tangent to the geodetic meridian of P.) Positive 
in the direction of increasing <p . See Fig^es 2.1 and 2.2. 

iVt it Z) Local Ellipsoidal Frame 

Origin: At the point P ( X, , u ). The above definition for the topo- 

centric frame applies here also. 
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C axis: Normal to the confocal ellipsoid of semiiniiior axis u , which 

passes through P. Positive in the outward direction. 

17 axis: Normal to ^ and the ellipsoidal meridian plane of P (i. e. , 

tangent to the ellipsoidal parallel of P. ) Positive in the 
direction of increasing X . 

4 axis: Normal to r) and the confocal hyperbolo.d passing through P 

(i.e., tangent to the confocal ellipsoid at P. ) Positive in the 
direction of increasing j 3 . 

(r)*, i*, C*) s Local Astronomic Frame 
Origin: At point P. 

C* axis: Normal through P to the geop of P (i. e. tangent at P to the 

plumb line passing through P.) Positive outwards, 
rj* axis: Normal to C* and to the mean astr momlcal meridian <'i p 

(positive in the direction of increasing astronomic ’ou^itude.) 
axis: Normal io C* :'nd r)* forming a right-handed system. I’ositive 

in the direction of increas'.ag astronomic longitude. 

In the ideal case of parallelism between the (x, y, z) and (u, v, w) 
systems, the transformation b /een the astronomic and geodetic 
coordinates ' done through the deflection of the vertical components 
17' and 4'- 
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3. 


CURVILINEAR GEODETIC COORDINATES 


3.1 General Comments 


Consider three families of surfaces represented by the parametric equations: 

X = X(u, V, w) <p = <p(u, V, w) h = h(u, v, w) (3.1-1) 

where (u, v, w) are Cartesian coordinates. This is really a transformation between 
points in the (u, v, w) Euclidean space lying in a certain domain and points in a 
certain domain in the (X, <p, h) space, generally a "non-flat" space. These domains 
naturally will exclude all singular points of the transformation. 

Assuming now that X, <p, h are variable param^ers, for each constant value 
of the parameters the family of surfaces will define three "coordinate surfaces" inter- 
secting in "coordinate lines or curves", hi general, one surface of each family passes 
throu^ a chosen point and a neighboring point will be determined by nei^boring values 
of the parameters, thus dividing the space into elementary cells which in general are 
not rectangular parallelepipedons. ff to each value of (u, v, w) corresponds a unique 
value of (X, (O, h), then any point P is uniquely determined by the three surfaces 
throi^h the point. 

The quantities X, <p, h are called the "curvilinear coordinates" of the point P. 

The most convenient system of curvilinear coordinates for geodetic fq)plica- 
tions are determined by families of surfaces which intersect each other everywhere 
at right angles, hi such a case we have a "triply -orthogonal" family of svmfaces or 
an "orthogonal curvilinear system". 

Assuming that (3.1-1) represents any set of orthogonal curvilinear coordi- 
nates, in that which follows, the general theory is going to be particularized, first 
to the very well-known set of geodetic coordinates and later to some other curvilin- 
ear orthogonal systems used frequently in geodesy and geophysics. All the matrix 
relationships derived, even though deduced for a particular curvilinear geodetic 
system, may be applied to :my set of orthogonal curvilinear coordinates. 
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3.2 The Local Base ( ei , e?, e O 


The coordinate transformation b^ween the curvilinear geodetic coordinates 

and the Cartesian coordinates may be expressed symbolically by 

(X, <p, h) (u, V, w) 

U,f ) 


and is defined by the well-known matrix relation 


u 


(N +h)co8<pcosX 

V 


(N +h)cos<psinX 

w 


(N(l - e^ ) + h] sin<p 


(3.2-1) 


where N, the principal radius of curvature in the prime vertical plane is given by 


N = 


(1 - e^siri^ (p) 


TTZ 


(3.2-2) 


and 


? _ 


= 2f - f 


exactly. 


Relation (3.2-1) can also be expressed in general by the usual parametric form: 
u=u(X,(p, h) v = v(X,(D, h) w=w(X, <p, h) (3.2-3) 

In order to have a coordinate system of practical value, the following conditions will 
be satisfied everywhere except at isolated singular |K>ints (e.g. , the polos): 

a) Each point (u, v, w) has a unique set of curvilinear coordinates: 
that is, there is a one-to-one correspondence between the (X, <p, h) 
and (u, V, w) coordinates. Therefore, the Jacobian determinant of 
truiisformation (3.2-3) is not zero. 

b) Equation (3.2-3) can be solved for X, <p, h giving the inverse transfor- 
mation: 


X = X(u, V, w) (p = (P(u, V, w) h = h(u, v, w) (3.2-4) 

This cannot be obtained by an explicit simple closed expression, but can 
be implemented through iteration [Heiskanen and Moritz, 1967], [Rapp, 
1975], [Bartelme and Meissl, 1975] or directly (Paul, 1973] and [Benning, 
1974]. 
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d) In addition to these conditions, the tangents at a point P to the X, <p, h 
"coordinate lines" through this point are perpendicular, so that the 
curvilinear system is orthogonal. 


Equation (3.2-3) in vector notation may be written as: 


where 


r = r(X, (p, h) 


r = ui + V j + wk 


(3.2-5) 

(3.2-6) 


and i, j, k are the unit vectors along the u, v, w axes. 


The tangent vectors to the (curvilinear) coordinate lines at P are defined by: 


hr 

Bu 

dX 

:♦ dv :» dw ;♦ 
‘ ' dX^ ' dX^ 



(3.2-7a) 

hip 

bu 

d(p 

:» dv hw ;* 

1 ' — J ^ — k 

hip hip 



(3.2-7b) 

dr 

dh 

du 

bh 

:♦ dv ? dw 



(3.2-7C) 

From (3.2-1) one can obtain 






^ = -(N h)cos<psinX 

o A 

b u 

b(p 

-(M t h) simp cos X 

du 

dh 

COS0CO8 X 


9v . 

r-T - (N + h) cosipcos X 
6 A 

bv 

b0 

-(M t h)sin<psinX 

dv 

dh 

cos (p sin X 

(3.2-8) 


bw 

b0 

(M + h)cos(p 

dw 

dh 

sincp 



where M, the principal radius of curvature in the meridian plane is 


Usii^ (3.2-7) and (3.2-8) 
orthogonal, that is. 


M 

it can 


a(l - e^) 

(1 - e sin <p) 

be shown that the 


(3.2-9) 

Cl (i 1, 2, 3) vectors arc 


Cl • c, 0 Vi / J 

Computations involving curvilinear coordinates are 


(3.2-10) 

greatly simplified if the 
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coordinate curves and the vectors c, are orthogonal, as in this case; otherwise the 
introduction of tensors will be required. 

From (3.2-7) and (3.2-8) it can be seen that the Cj (i ^ 1, 2, 3) vectors are not 
unit vectors; thus it will be practical to replace the Cj by unit vectors e^ (i = 1, 2, 3) 
having the same directions. 

Defining, 


ei 


^2 




1 d r 
hi 

1 ^ 

ha <10 

1 ^ 

ha ^h 


(3.2-lla) 

(3.2-llb) 


(3.2-lIc) 


where hj represents the corresponding modulus of the vectors Cj, or: 


hi 


r^rl [/duV /•^v\‘ 

dh I ~ ?)h/ ^\^^h ) ' \dh ) J 


^ I 
dX 
-> 

50 


(N t h)cos0 


(M ' h) 


(3.2-12a) 

(3.2-12b) 


(3.2-12C) 


Thus using expressions (3.2-7), (3.2-8) and (3.2-12) in (3.2-11) 

-> - * 

Cl - sinX i ' cos Xj 

Ca -sin0cosXi - sinosinXj ' cos 0 k 

Ca - cos0cosXi + cos0sinXj ^ sin 0 k 


(3.2-13a) 

(3.2-131)) 

(3.2-13C) 


ft can be observed that the vectors e, (i 1, 2, 3) arc of unit length and mutually 
orthogonal, i.e.. 




6i, 


n (i j) 

< 0 (i / j) 


(3.2-M) 


where 6ij in the above formula is referred to as the "Kroneeker delta". 

-> 

'I'hus the vectors ci form an orthonormal base in the Kuclidean spaje ic' just as 
do the vectors i, j, k. There is, however, one fundamental difference between the two 
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bases: whereas the i, j, k are fixed directions, the directions of et(i - 1>2, 3) in 
general will vaty from point to point because the coordinate lines are curved. This 
can be seen h'oi i (3. 2-13) where the ej vectors are functions of X and (p. A frame 
like (ei, e-j, which is parameter-dependent is called a "movit^ frame." The theory 
of the "msvin frames of reference" was greatly extended and generalized by [Cartan, 
1935] who referred to it in French as "rephre mobil. " For an introductory stud|y of 
the matter aud some bibliography, see [ Grafarend, 1975]. 

The rei tangular Cartesian reference frame (t}, 4* D whose axes have the same 
direction as tt : unit vectors et is said to be the "local coordinate reference 
frame' which ; s attached to the point P. The coordinates (or components) of at(y vector 
V in this local movii^ frame (T), 0 are termed the local coordinates of v. Thus, 


P = t?ei + |ea +Ce3 


(3.2-15) 


K is important to note that because the vectors et are functions of X and <p, in general 
the compon >nts of ^ are not ^ ^ ^ but 

ax ax’ ax’ ax 

av_ hr) -* ^ ^ f ^62 ^ . a eta o \c\ 


ax ^ax ^ax 


hV 

Clearly the same logic wi’ apply to — . This dependence of the vectors Ct on the geo- 

a<p 

detic coordinates will be implied always when a vector u is written in the form u(r), C* 0* 


? . 3 The Rotation Matrix R 


Denoting by [u v w] ' the 3x1 matrix whose elements are the coordinates 
— > _ 

of a free vector v in t..^ fixed Cartesian frame and by [77 i C] the column matrix 

whose elements : "e the coordinates of 1/ in the local system, the transformation of 

components of the free vector l> from a reference frame to the other is given by; 


T) 


u 


U 

V 

r 


w 

# * 
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or symbolically may be represented by the mapp'ng 


(u, V, w) 


-> (T|, 0 


where the rotation matrix R of the transformation can be deduced from Fipire 2.2 by 
simple geometric considerations as follows: 


R = Ri(90 - o) Ra(X 90) 


” sinX 

cosX 

0 

- sinpeos X 

- sinipsinX 

cosio 

cos (0 cos X 

cosipsinX 

sin0 


(3.3-2) 


H; must be pointed out here that in the following the interpretation of R will always be 
that of an orthogonal transformation from the geocentric to the local system. Neverthe- 
less, knowing that R is an orthogonal matrix (RR^ I > It' R^), the inverse trans- 
formation can also be written as 


(T). 0 


it' 

(u, V, w) 


ft should be noted that the rows of R are the components of the orthonormal 
vectors ei (i 1, 2, 3) given in (3.2-13). Therefore it follows immediately that 




— ^ 



* 





It 

j 



—V 

C3 


k 

L J 


which gives the transformation between the two orthonormal bases, ft is known that 
the rows of R represent as well the direction cosines of the vectors e^ (i 1, 2, 3), 
that is, the direction cosines of the normals with respect to the throe ooordimiLe lines 
or surfaces at P. 


3.4 Element of Arc and Orthogonality 

ft was proved already’ that the vectors e, (i 1, 2, 3) tangent to the coordinate 
lines are orthogonal. 

A different wav to see that the curvilinear geodetic coordinates are orthogonal 
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is computing the element of arc (element of distance) in these coordinates. From 
(3. 2-3) 


u = u(X, <p, h) V - v(X, <p, h) w = w(X, 0 , h) 
The total differentials of the functions u, v, w are 


. . . du , Su .. 

du = :^dX^ — d0^ — dh 
dX 3 h 

, dv j - dv , 9v 

dv :-r dX ^ — d 0 + ^ dh 

ax S0 ah 

. 3w aw j Sw .. 

dw — dX ^ — d 0 + ^ dh 

ax o0 ah 


or in matrix notation 


du 


dX 

dv 

- J 

d 0 

dw 


dh 

J 


where J is the Jacobian or functional matrix and may be expressed by: 


(3. 4-la) 
(3.4-lb) 


(3.4-lc) 


(3.4-2) 


j = w) 

a(X, 0 , h) 


{ ^r w ) 
(X, 0 , h) 


“ du 


du 

ax 

d0 

dh 

dv 

dv 

dv 

dX 

d0 

dh 

dw 

dw 

dw 

_ dX 

d0 

dh 


(3.4-3) 


Thus the Jacobian matrix of a coordinate transformation can be interpreted as the 
matrix of a certain linear change of coordinates, namely: 

(dX, d0, dh) — — > (du, dv, dw) 

As mentioned previously, the Jaoobi:m dotorni inant is not zero. Therefore 
I J I / 0 > J is nut singular > J‘ exists. 

The square of the lino element in the (u, v, w) system is 

ir> 



ds" = du^ +dv^+dw^ 


(3.4-4) 


Substituting above the values from (3.4-1) one can obtain the following equation 
ds^ = hx^dX'+ ha^d<p + ha dli + 2 h 4 "dXd<p + 2 hfe^ dcpdh + 2 he‘ dhdX 
where the values hi, hg, ha are given by (3.2-12) and 


du 

Su 

_L 

3v 

3 V 


3w 

3w 


dX 

do 

r 

dX 

30 


3X 

30 _ 


ii? 

du 



3v 


3w 

3w 


30 

dh 

r 

30 

3h 


30 

3h 


3u 

du 


3 V 

3 V 


cw 

3w 


3h 

dX 


3h 

3X 


3h 

3X . 



(3.4-5) 

(3.4-6a) 

(3.4-6b) 

(3.4-6C) 


Replacing the values presented in (3.2-8) in equation (3.4-6), it is easy to find that the 
necessary and sufficient condition for orthogonality is 


h4 = he = h« = 0 (3.4-7) 

Therefore the absence of the terms dXd 0 , d 0 dh, and dhdX in (3.4-5) is the 
evidence that the curvilinear coordinates (X, 0 , h) are orthogonal. The transforma- 
tion is called conformal when elements of arc in the neighborhood of a point In the 
(u, V, w) system are proportional to the elements of arc in the neighborhood of the 
corresponding point in the (X, 0 , h) curvilinear system. That is when 

ds"^ = du"" + dv^ +dw^ = k^(dX" +d</f +dh") (3.4-8) 


Thus 'M formality requires 

hi - hg = ha (3,4-9) 

The above is in agreement with the fact that conformality implies orthogonality but not 
viceversa. 

ft can be observed that equation (3.4-5) gives the linear element ds In three- 
dimensional space. Clccirly when limited to surface transformations the so-called 
Gaussian fundamental quantities will be present. For example, in the case of an 
eliii.soid the following identities are established: 

K hi G = hg F = I 14 (3.4-10) 
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When F = 0 the condition for orthogonality exists, and if simultaneously E = G 
the transformation is conformal. 

Equation (3.4-5) can easily be written in matrix notation as follows: 


dS^ = [du dv dwl 


du 

dv 

dw 


But recalling equation (3.4-2) 


( 3 . 4 - 11 ) 


ds'" = (dX d<p dhlJ^J 


dX 

d<p 

dhj 


(3.4-12) 


Thus the transformation between two sets of coordinates in E^ will be orthogonal if 
the matrix product of the Jacobian transpose by the Jacobian is a diagonal matrix. 


3.5 The "Metric Matrix” H 


The differentials du, dv, dw may also be represented as a free vector. They 
behave under transformation of coordinates as do free vector components; thus from 
equation (3.3-1) 


dt) 


du 


= u 

dv 

1 


dw 


Using (3.4-2) it is possible to express the relationship between the differential 
arc-length of the curvilinear coordinates along the "coordinate lines" and their pro- 
jection on the (t), ^ , C) system as 
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(3.5-2) 




dX 

d4 

= RJ 

d<p 

dC 


dh 


v'here the value of the matrix product ItJ according to (3.3-2), (3.2-8) and (3.1-3) is 



(N +h)cos<p 

0 

0 


hx 

0 

0 

RJ = 

0 

(M +h) 

0 


0 


0 


0 

0 

1 


0 

0 

ha 


- 


_ 


- 


_ 


Thus symbolically one may represent 

(dX, d<p, dh) (dT7, d4, dO 


where the matrix of the transformation H will be called the "metric matrix" by analogy 
with tensor terminology. 

From (3.5-3) the following basic relation can be written; 

IIJ = H (3.5-4) 

and using the orthogonality property of U, 

J = H'H (3.6-6) 

Then 

yj = 

and H being diagonal the following results: 

J^J = U' (3.5-6) 

as can easily be proved by simple multiplication of matrices. 

Comparing (3.5-6) with (3.4-12) it may be deduced that the diagonality of the 
metric matrix H is a consequence oi tlic orthogonality of the curvilinear coordinate 
system. 
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3.6 Jacobian Determinant and fts Applications 


Equation (3.5-6) provides a simple way to obtain the value of the Jacobian 
determinant. Taking determinants in (3.5-6) 

I JTJ I = I hM => I JM I J 1 = I H P but I jM = I J 1 

Thus 

I J P = I H P 

and finally 

I J I = 1 H 1 = hi haha (3.6-1) 

where hi (i = 1, 2, 3) are given by equation (3.2-12). 

B; is important to mention here that while the determinants of J and H are always 
equal, the matrix J is equal to the matrix H only when R = I, as can be seen &om (3.5-4). 
This will be equivalent to making the frame (u, v, w) parallel to the local (r), C) 
frame through the pertinent rotations. 

The functional (or Jacobian) determinant in the case of geodetic coordinates is 

I J 1 = (M +h) (N +h)cos<p (3.6-2) 

Thus, aside from points where cos0 = 0 =^<p = ±|- ^1J| = 0, the transforma- 

tion is locally one-to-one, implying that any point on the w (polar) axis is a singular 
point in this specific transformation. 

Equation (3.6-1) is also very convenient for the computation of elements of area 
along the different coordinate surfaces and the element of volume between them. 


dAi = hihadXd(p (3.6-3a) 

dAe = hihadXdh (3.6-3b) 

dAa = hah3d<pdh (3.n-:ic) 

dV = hihphidXd<pdh (3.6-3d) 

For the paxlicular case of geodetic coordinates on the reference ellipsoid 

dA - MNcosodXdca (3.6-4a) 

dV = (M +h)(N +-h)co8todXd(pdh (3.6-4b) 
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ft should be noted that the Jacobian may be either positive or negative, the 
difference of sign being of the same nature as the consideration of an area or volume 
and their reflection. Due to the fact that the elements of area or volume are con- 
sidered positive when the incremenia of the variables are positive, the absolute value 
of the Jacobian will always be taken. 


3.7 Analytic Hxpressions for the Inverse of the Jacobian Matrix 


From (3.5-5) it is evident that 

J^ = hMi (3.7-1) 

which provides an analytical way of obtaining the inverse of the Jacobian matrix. Using 
(3.5-3) and (3.3-2) one obtains 


J^ 




* ^ 


^ - 



d\X 

dv 

dw 

a(X, <p, h) 

. (X, (p, h) 

5(p 


dtp 

d(u, V, w) 

(U, V, w) 

du 

d V 

dw 









d V 

dw 


sinX 


cos X 


(N +h)cos<D 

siniQcosX 
M +h 

cos (p cos X 


(N + h) cos<0 

sin<psinX 
M +h 

co8<p sin X 


0 

C 08 <p 

M +h 
sin(p 


(3.7-2) 


The classical way to obtain the elements of J ^ Is to solve the following nine 
equations: 


du 


+ ^ 


+ 

1!} 

dX 

du 

d<p 

du 

dh 

du 

dV 

dX 





dX 

dv 

dtp 

dv 

dh 

dv 


(3.7-3a) 

(3.7-3b) 
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5w 



aw 


4. 

aw 


dX 

aw 


a<p aw 

T 

dh 

aw 

5u 



au 



au 


dX 

av 

1 

acp 

a V 


dh 

av 

5u 


+ 



+ 

a u 

i!i 

ax 

aw 


a(/5 

aw 


ah 

3w 

av 



av 


4. 

av 


ax 

au 


a<p 

au 


ah 

au 

av 



a v 


4. 

av 


ax 

aw 

T 

a0 

aw 


ah 

aw 

aw 


4 * 

aw 

a<p 

4. 

aw 

Hi 

ax 

au 

r 

b((> 

dU 


ah 

au 

aw 



aw 

bip 

4. 

aw 


ax 

av 

1 

ao 

av 

r 

ah 

a V 


which reduce in matrix notation to 

JJ‘ = I 


(3.7-Sc) 


(3.7-3d) 


(3.7-3e) 


(3.7-3f) 


(3.7-3g' 


(3.7-3h) 

(3.7-31) 


(3.7-4) 


or 

(jj^)^ = I =» (j')^ / = I 


(3.7-5) 


but from (3.7-1) == , thus substituting this above, the following two equalities 

can be written: 

R^H'J^ = I (3.7-6) 

and 

JH‘U = I (3.7-7) 


Premultiplying both sides of (3.7-6) by j\ which is equal to H’^R, 

(H"R) (R'h‘)J^ = J‘ or 

= (H") ^3.7-8) 

where clearly (H*)^ = (H‘^) ‘ due to the diagonality of H. This Is another way of 
computing independent of R as a function of J and H. Substituting (3. 7-H) in 
(3.7-4) the following equality is established: 
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J(H"^) = I 

Equation (3.7-8) gives 

J = [(H“)V]‘ = (J^)‘H“ = (J^)’H“ 

Therefore 

J = (3.7-9) 

3.8 The Matrix R as a function of H and J 

Finally it is possible to write the expressions for the rotation matrix 11 as a 
function of the md:ric and Jacobian matrices. 

From (3.5-4) 

U = IIJ' (3.8-1) 

Substituting (3.7-8) above, the following is derived: 

II = (3.8-2) 

and therefore the following matrix equality holds: 

HJ'- H'J^ (3.8-3) 

Knowing, as was mentioned abi ve (section 3.3) that the elements of the rows of 11 
represent the direction cosines of the normals to the family of surfaces X, <p, h 
that pass through a point P, equations (3. -li and (3.8-2) will provide genera' *‘orm- 
ulas for computii^ the nine direction cosines of these normals in two different ways. 

For example, the direction cosines of the normal to the reference ellipsoid 
‘it a point (X, <p, h) may be given by 

1 du 8h 

cos (u, h) - ~ rr" ~ ha — = coscpcosX 

ha o n o u 

/ i-v 1 dV . dh . . 

cos (v, h) ^ ^ ^ "" cososmX 

... 1 aw .ah 

cos(w,h) = - _ . h, _ - ta® 

This approach will be very helpful for other curvilinerjr orthogonal systems where the 
matrix U cannot be obtained by simple geometric considerations, as in the case of 
curvilinear ellipsoidal coordinates (see Section 5. 2). 



4. 


DIFFERENTIAL CHANGES BETWEEN CARTESIAN AND 
CURVILINEAR GEODETIC COORDINATES 


4.1 Basic Equations 

From (3.5-2) the following basic relation can be written: 


dT)‘ 


’ dX' 

dC 

= H 

d(p 


1 

1 

dh 


(4.1-1) 


which substituted in (3.5-1) gives the fundamental differential relations b^ween the 
Cartesian (u, v, w) and curvilinear (X, <Pt h) coordinates. 


dX 


du 

d(p 

= II 'R 

dv 

dh 


dw 


(4.1-2) 


This is usually written for the case of geodetic coordinates as 



dX 


du 

H 

dtp 

= R 

dv 


dh 


dw 

- 


(4.1-3) 


or 


(N +h)cos(pdX 


du 

(M + h) do 

= II 

dv 

dh 

(•.n 

dw 

L. 


(4.1-4) 


where a particulai ellipsoid is implied in the computation of N and M. 

The above formula expresses the basic matrix equation relating the differen- 
tial changes in the geodetic coordinates (dX, d(/J, dh) of a point P referenced to a given 
ellipsoid (a, f), d<'.e to differential changes in the geodetic Cartesian coordinates of 
the point. 
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Considering the orthogonality of R 


du 


- 

(N + h)cos(pdX 

dv 

= R^ 

(M + h) d<p 

dw 


dh 


ft is possible to show that this equation is equivalent to (3.4-2), namely. 


du 



’ dX' 

dv 

1) 

V, w) 
<p, h) 

d(p 

dw 



dh 


Clearly, if the point P(X, <o, h) is on the surface of the ellipsoid 


du 


Ncos( 0 dX 

dv 

= R' 

Md(p 

dw 


0 


(4.1-5) 


(4.1-6) 


4.2 Differential Changes in (X. <o, h) Due to Siifts, Rotations and Scaling of the 
(u. V. w) Cartesian System 

As an illustration of the above theory, one can assume, for example, that it is 
desired tu obtain the differential changes in the geodetic curvilinear coordinates 
(X, <p, h) due to differential shift , rotation and scale changes of the Cartesian 
system. Then from (4.1-4) the total contribution may be expressed as: 


(N +h)cos<pdX 



du 


du 


du 

j 

(M + h) do 

= R 


dv 

+ 

dv 

+ 

dv 


dh j 

1 


d w 

shift 

dw 

rotftt ion 

dw 

•QftlQ ' 


Each Individual differential contribution will be studied separately in the following 
paragraphs. 
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4.2.1 Changes Due to Differential Shift of the Origin 


ff the geodetic sy^em (u, v, w) is shifted by the amounts 6 u, 6 v, 6w> then 
obviously one obtains 


du 


5u 

dv 

= 

6v 

dw 

shift 

6w 


Assiuning the systems (u, v, w) and (x, y, z) to be parallel, the signs of the shift 
components may be given by one of the following conventions (see also Fig. 4.1): 




’du ■ 


dv 


dw 



(Geographic i^stem) - (Geodetic i^stem) 
(Final " ) - (Initial " ) 

(New " ) - (Old " ) 

(Fixed ” ) - (Moving '* ) 


(4.2-3) 


Thus for example the transformation of coordinates 

Geodetic > Geographic will be 

(u, V, w) (X, y, z) 

Geographic = Geodetic + (Geographic ~ Geodetic) 


and consequently 


X 


u 


6u 

y 

= 

V 


6v 

z 


w 


6w 


(4.2-4) 


In general the above sign rules will be observed on the following pages, if not specified 
otherwise. 


4.2.2 Changes Due to notation 

As is known, it is possible to relate the coordinates of the two Cartesian 
systems having the same origin by the equation 
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(4.2-5) 


- 

X 


U 

y 

= B 

V 

z 


w 


where the rotation matrix R can be written as follows 


R = K^(i(i)Ri(€)R?(aJ) = 


cosiCcos^ - sino)smesin(^ 
- sinc*.'cop t 

cos cosing sinoisin ccos(() 


sin cocos ^ + cos cosine sin ^ 
cos cocos e 

sincosind) - cos cosin e cos 


- cose sin 
sin ^ 

cos € cos ^ 


(4.2-6) 


An introductory section on rotational matrix algebra may be consulted in [Goldstein, 19SU]. 

In order to keep the sign convention established in the previous section, the dif- 
ferential chaises in the coordinates (u, v, w) due to the rotations co, (, t{) are given by 


du 


X 


u 


- 

u 


u 


u 

dv 

= 

y 

- 

V 

= R 

V 

- 

V 

= [R- I] 

V 

dw 

rota t too 

z 


w 


w 


w 


w 


(4.2-7) 


Assuming now differentially small rotations 6co, 6c, 6(|) (see Fig. 4.2), it 

follows 



sin CO 

3 

<o 





sin c 

« 6c 



(4.2-8) 


sint(j 

ss 6(^ 





cos CO 

= cos € 

= CO8 0 

^ 1 


and neglecting second-order terms, the rotation matrix becomes 



1 

6co 

- 6il) 



U 

- 6co 

1 

6 c 

= 

(1.2-9) 


6ip 

- 6c 

1 
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and by means of (4.2-7) 


du 


0 

6co 

- 6(() 


u 


“ — 
u 


dv 

= 

- 6w 

0 

6 ( 


V 

= 6R 

V 

(4.2-10) 

dw 


6(|j 

- 6? 

0 


w 


w 



B; can be observed that the matrix 6 R is a skew-symmetric (or antisymmetric) matrix 


6R+(6R)" = 0 

That is, if 6rij (i, j = 1, 2, 3) are the elements of 6 R 

6rij =0 V i = j 
= - 6rij V i j 


(4.2-11) 


(4.2-12) 


This represents an important property for all differential rotation matrices. 
Notice that the R& matrix given by (4.2-9) is orthogonal only up to first-order terms. 
From (4.2-9) and (4.2-10) 


Ra = I +6R 


(4.2-13) 


The condition for Rato be orthogonal is RaRa’^ = I . but 

RalV = (1 + 6R) (1+ 6U)' = (I ^ 61^ (I- 6R) = I- (6R) ' ^ I 


4.2.3 Changes Due to Differential Scale Changes 

Ihe changes in u, v, w due to a differcnti:il scale change 6 1. are 

(4.2-14) 

where 6L according to the sign convention mentioned in (4.2.1) is, 

6L (Geographic scale) - (Geodetic scale) (4.2-15) 

Kquation (4.2-14) can also be written as 


du 


u6L 


u 

dv 

= 

v6 L 

- 6L 

V 

dw 

5l 

w6 L 


w 
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(4.2-16) 


X 


u 


u 

y 

- 

V 

= 6L 

V 

z 


w 


w 


From (4.2-16) it is obvious that 


X 

y 

= (1+6L) 

u 

V 

z 

— ^ 


w 


(4.2-17) 


which implies that the 6 L change in scale may be computed directly from the coordi- 
nates, but 6L can also be obtained from the chord distances d^ and dy in the respective 
systems. 


d/ = (x - x')" + (y - y')^ + (z - z')^ = (x - x' y - y* z - z*] 


- = 


X - x’ 

y -y* 

z - z' 


= (1 + 6L)^ (u - u* V - v' w - w'] 


u - u’ 


V - V* 


w - w* 


= (l+6L)‘’dy (4.2-18) 


Thus 


d* = a +6L)du => 6L = 

dy 


(4.2-19) 


Precautions should be taken, however, when in a least square adjustment, scale d^er- 
mination through chord distances is intended. Only an independent set of chords 
should be used [Leick and van Gelder, 1975). 

4 . 2.4 Final Equation 


Thus finally substituting in (4.2-1) the computed effects in the coordinates 
u, V, w due to shift , rotation and scale changes, the following equation is obtained: 
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(N + h)cos<pdX 



- 

6u 


0 

6o) 

-60 


u 


u 

) 

(M + h) d<p 

= R 


6v 


- 6co 

0 

6c 


V 

+ 6L 

V 


dh 

_ 

1 


6w 


60 

-6c 

0 


w 


w 

i 


(4.2-20) 


which gives in matrix notation the changes in geod^io coordinates <A» c 3 , h) due to shifts 
(6u, 6v, 6w), rotations 6u), 6f, 6^ and scale 6L 

Clearly to express dX, d(p, dh only in ftinction of geodetic coordinates, equation 
(3.2-1) will be substituted in (4.2-20). 


4.3 Similarity Transformations 


ft is proper to point out here that if one considers only changes in the Cartesian 
coordinates due to translations, rotations and scale change, after substituting (4.1-4) in 
the left side of (4.2-20) the following results: 



du 

1 


6u 

R 

dv 

= R 


6 V 


dw 


i 

6w 


0 6co -60 

- 6<j 0 0 6 f 

60 - 6 € 0 




- 

. 

u 


u 


V 

+ 6L 

V 


w 


w 

1 


But according to the notation used in this report. 


du 


X 


U 

dv 


y 

- 

V 

dw 


z 

f tnftl 

w 

L J 


(4.3-2) 

Jlnlttal 

Thus, finally after substitution of the above in (4.3-1) and omitting the rotation matrix it 


X 


u 


6u 


0 

6u) 

- 60 


u 


u 

y 


V 

+ 

6v 


- 6C0 

0 

6c 


V 

+ 6L 

V 

z 


w 


6w 


60 

_ 

- 6c 

0 


w 


w 


(4.3-3) 


But one can also write 


6 u = Ax 

6v = A = i S®odetIc systena with 

^ ^ ~ I I espect to the geographic system 

6w = A z 

5L = A ^ Geographic scale - Geodetic scale 


(4.3-4a) 

(4.3-4b) 


Then one obtains the usual seven-parameter transformation between the coordinates 
uf two Cartesian frames. 



Some authors use the notation 1 + A = A . Thus, recalling (4.3-7) 


(4.3-5) 


(4.3-6) 


(4.3-7) 



which is a form of the more general similarity transformation: 


(4.3-8) 


X = T + ARU 


(4.3-9) 


Comparing (4.3-8) with (4.3-9) it is obvious that besides the assumption of small rota- 
tions R « R* the products of A by the rotations were neglected in the differential 
transformations of the formula (4.3-8). For a complete development of the similarity 
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transformations as used in geodesy* consult the discussion in [Leick and van Gelder, 
1975]. 

Any one of the above type transformation models is generally referred to in 
geodetic literature as the "Bursa model" after [Bimsa, 1966]. They have been popu- 
larized as mathematical models of the type F(X, L) = 0 in the least square solution 
for computit^ the seven parameters of the similarity transformation between world 
systems. 

Nevertheless, it is unclear why several authors, among them [Badekas, 1969] 
and [Krakiwslqr and Thompson, 1974] credited {Molodenskii et al. , 1962] with a dif- 
ferent model where the rotations and scale erpansion are about some particular point 
(Uo» Vo, Wo) other than the origin of the Cartesian geodetic system. 

hi the following sections the equations given by [Molodenskii et al. , 1962] will 
be presented adopting the general criterion of this work, h will be shown that they 
are not different from a similarity transformation of the form (4.3-5) except in 
the way the scale is applied. In th connection the similarities of eliminating the 
variations of scale through changes in the semi major axis of the reference ellipsoid 
will be explained. 

4.4 Differential Transformations According to [Molodenskii et al. , 19621 

Before fully developing the formulas given by [Molodenskii et al. , 1962], 
thf effect of differential changes in the Cartesian geodetic coordinates due 
to changes in the size (a) and flattening (f) of the ellipsoid will be 
treated. 


4.4.1 Changes Due to Variations 6a and 6f 


If the original reference ellipsoid parameters are changed, their effect on the 
(u, V, w) coordinates can be expressed In matrix notation as follows: 


where 


6a a 

6f = OM 

and from (3.2-1) it is possible to compute (Uapp, 1975a]: 


where 


f)U _ C03<p COB X 

3a " W 

3v _ cos<psinX 
da “ W 

^ _ (1 - e^) sin0 
3a “ W 


= 1 - e“sin“<p 


3u _ a(l - 0 8in“<pco3<peo8 X 
3f ^ 

3v a(l - 0 ain^VcoscpslnX 
3f W’ 


— = (Msin“<p- 2N) (1 - f)sin<o 

O I 


m=2S^ 


N = — 
W 


4.4.2 General Equations 

Literally following Molodenskii et al., [1962], their equation (1.3.2) using 
the ciiffirential matrix approach and the notation of this paper may be written as 

du 6uo 0 6co - 6([) u - Uo 

dv = 6 Vo + - 6co 0 6 c V - Vo 


w- Wo 


+ R" 


(N +h)cos<pdX 
(M + h)d0 
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Fig. 4.3 Datum Origin Coordinates in the Geographic 
and Geodetic Sjrstems 
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Clearly in the above equation differential changes (dX, dtp, dh) in the geodetic coordi- 
nates at any particular point are considered, in addition to the changes 6 a and 6f. 

Before goin^ further, it is important to understand v’hat [6ua 6vo 6wJ^ is. 
This is not very clear in the original text and probably originated the confusion when 
equation (4.4-5) was used as a model. (At this moment it is proper to mention that 
in the often quoted work by Molodenskii et al. , [1962] a local system different from 
the one used in this paper is assumed. Appendix A contains the relationship between 
the matrix notation used *'ere and the equations in the original English translation from 
the Russian). 

Followin, [Molodenskii et al. , 1962] the shifts between the origins of the geo- 
graphic and geodetic systems Vecome (see Fig. 4.3) 


du 


1 

o 

L 


Uo 


Uo 



Uo 


dUo ^ 


dv 

= 

Vo 



= 

<1 

o 

- 


Vo 


dVo 


dw 

1. 

X 

1 

o 

_ 1 

X 

_ '^0 . 


o 

1 

X 1 

> i 

Wo 

u 

dWo 

rot / 



Uo 



Uq 


0 

6co 

- 60 


Uo 


= 

Vo 

- ' 


Vo 


- 6co 

0 

6c 


Vo 



. 

X i 


Lwo_ 

U 

60 

- 6c 

0 



u 


where for clarity , small subindices are used to represent the Cartesian system to 
which the components of the column vectors are referred. Consequently, it is possible 
to write the following relation equivalent to equation (1.3.4) in [Molodenskii et al. , 1962] 
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and therefore from (4.4-6) and (4.4-7) one may conclude 
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(4.4-7) 
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(4.1-8) 



Probably the main reason why in [Molodenskii et al. , 1962] the value of (4.4-8) is not 
given explicitly is that what th^r called "progressive translations of the ellipsoid" are 
the difference between the two sets of coordinates of a particular point (in this case 
the origin of the datum) which are in different Cartesian systems (x, y, z and u, v, w). 
This difference, although rigorously correct, is somewhat difficult to visualize. 
Nevertheless it is perfectly clear that equation (4.4-8) does not represent the shifts 
between the origins of the (x, y, z) and (u, v, w) systems, as interpreted by some 
authors when the Molodenslui equations were used as a model. It is obvious that the 
second term on the right hand side of (4.4-7), equivalent .o a rotation about (Ug, Vg, Wg) 
of a system parallel to (u, v, w), is neglected when the vector [6ug Sv^ 6wg|^ is 
thought to be the shifts between the geographic and geodetic Cartesian systems. Sub- 
stituting the value of (6Ug 6Vg 6 Wg]’^ from (4.4-7) in (4.4-5) It may be shown that 
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(4.4-9) 


and finally with the assumption (d X = d (0 = dh = 0), it follows that 
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(4.4-10) 



which is exactly the similarity transformation given by (4.3-5) with the notatioa 


du = Ax 
dv = Ay 
dw = Az 



shifts between the origins of the two coordinate systems 
with the sign convention of (4.3-4a) 


and the replacement of the effect of a scale change A by the change in the semi- 
major axis of the ellipsoid, hi the above equation it may le assumed 6 f = 0. 

Equation (4.4-10) is not given in the above form in [Molodenskii et al. , 1962] and was 
derived here only with the intention of showing that the model obtained is strictly 
similarity transformation of the "Bursa type" without anything special 
introduced besides the scaling variation mentioned above. Later the difference 
between these two scaling approaches will be fully explained. 

What Molodenskii presented as expression (1.3.3) is nothing else but equation 
(4.4-9) premultiplied on both sides by the rotation matrix It, namely 
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6 f 
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dw dw 
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4.5 Comparison of Scaling Methods by Means of 6 L or 6 a 


Before giving the individual formulation for each method, it must be understood 
that a change 6L is always applied to the unit length of the Cartesian coordinate system 
involved. Therefore it may be considered as scaling the space or a change in 
its metric. B; can be applied either to the geodetic or to the geographic system, 
depending on the adopted sign convention. 
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Nevertheless, a change 6 a is always applied to the geodetic system and prac- 
tically represents a "network scale. " 

4.5.1 Chaises dX, do, dh Due to a Change 6L in the Scale 

These changes may be obtained very easily by making use of the relations 
(4.1-4) and (4.2-14), namely 


(N + h) cosod X 


u6L 


(N +h)cosOcosX6L 

(M <- h) do 

= R 

v6L 

= R 

(N + h)cos0sinX6L 

dh 

6t 

w6 L 


(N(l - e“) +h]sino6L 


Finally, it can be proved that (see Appendix B) 


dXftt 

= 0 

(4.5-2a) 

d05u 

Ne^sinocoso 

STTh — 

(4.5-2b) 

dh5L 

- (aW +h)6L 

(4.5-2C) 


Consequently, as expected, for a reference rotational ellipsoid, there is not 
any influence in the geodetic longitude due to a change of scale in the length unit of 
the Cartesian geodetic system. 


4.5.2 Changes dX, do, dh Due to a Change 6a in the Semimajor Axis 
of the Ellipsoid 


The changes ^du, dv, dw) in the Cartesian coordinates in function of changes 
of the geodetic coordinates (dX, do, dh) and differential changes in ellipsoidal 
parameters 6 a, 6 f may be expressed by 




du du 



du 


da d f 
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(M ■<-h)cos<pdX 

dv 
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dv dv 

da df 



+ R^ 

(N + h) d<p 

dw 


dw dw 


0 I 


dh 


da df _ 




Thus, assuming that the point (u, v, w) remains fixed in space, the effect of 
differential changes Oa, 6f on (X, <0, h) is: 


(M + h)cos(pdX 

(N +h)d<p I = - R 

dh 


Therefore, recalling (4.4-1) and (4.4-3) 


du ^ 
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dv dv 
da d f 

dw ^ 
da df 
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6f 


(4.5-4) 


(M +h)cos(pdX 


du 

(N +h)d(p 

= - R 

dv 

dh 

5a 

dw 


R 


cos<pcosX . 

— «a 

co8<p8inX . 
W 


(4.5-5> 


The above matrix operations yield (see Appendix B) 


^’Xe. 

d06, 

dhfi* 


0 

Ne^8ln<pco8<p . 
(M+h)a 


- W6a 


(4.5-6a) 

(4.5-Gb) 

(4.5-6c) 


As in the previous case, changing the semimajor axis of the reference ellipsoid 
does not affect the geodetic longitude. Nevurtheluss, there are some differences 
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between the two methods of scaling, as will be explained in the next section. 

4.5.3 Comparison of Scaling Methods 

Assume that it is desired to absorb the metric scale change in the coordinate 
system by a change 6 a in the semimajor axis of the ellipsoid. A positive change 
6 L will not alter the ellipsoid; nevertheless, the new value of the semimajor axis 
(using the new yardstick due to a positive 6L) will fulfill the inequality, 

a < a 014 \/6L > 0 

Thus, 6ai_ ~ a ~ a oia ^ d, 

therefore, 6L > 0 =-> 6a^ <0, 

and finally, 6 at = -a6L (4.5-7) 

Therefore, substituting the above forn ula in (4.5-4b) and (4.5-4c), 

. N. v>3<p - 

d06«^ 6L = doftt (4.5-S) 

and 

dhft,^ = aW6L dhot = (aW ^h)6L (4.5-9) 

Thus when the semimajor axis of the ellipsoid changes by the amount 6 a,. 

due only to scale variation, namely (4.5-7), the differential changes in the latitude 

are equal to the differential changes produced by a scaling of the coordinate system. 

However, the differential changes in the heights are different, as ca i be seen from 

(4.5-9). This is in accordance with the remarks made by Hotine (1969, p. 264]: 

"Most of the systematic error in scale of a network could be 
eliminated by altering the size of the base spheroid in the geo- 
detic coordinate system. . . However, this procedure would 
vitiate the height dimension and would result in some inaccuracy 
even in a two-dimensional adjustment which ignores geodetic 
heights, especially if the network covers a considerable area." 

The difference between the two methods of scaling is 

Ah = dhftL-dha*^ = h6L (4.5-10) 
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which means that for points on the ellipsoid (h - 0) the two scaling methods are 
identical. 

ft can be observed that in any case the value of ^ h may be ne^ected, as can 
be shown by a simple example. Taking the values of 6L and h excessively, e.g. 

6L = 6 X 10~° and h = 5 km 

Ah = 30xl0‘*W. = 30 mm. 

Thus for all practical cases, the method of scaling followed by Molodenskii 
is not different from scaling the system throu^ changes in the unit length of the 
Cartesian axes. 

However, rigorously speaking, in general it is possible to assume changes 
in the semimajor axis of the ellipsoid 6a and at the same time, scale changes 6L. 

This, in fact, signifies that one can change the metric scale 6L of the Cartesian 
system without changing the size of the ellipsoid, although every length 
measured with the new scale unit will be different. This is thus equivalent to chain- 
ing the metric of the space, that is, its unit of length. 

On the other hand, it is also completely valid to assume a change 6a in the 
size of the ellipsoid independent of any scale change. This may be considered as a 
"network scale change" but clearly the imits along the axis of the Cartesian coordinate 
system will not undergo any variation. That is, the unit of length (i.e. the scale of 
the three-dimensional space) remains the same. 

4.6 Effects of Other Differential Changes on the Geodetic Curvilinear 

Coordinates 


The theory given in the previous sections is general and may be used in any 
case desired. The basic equation is (4.1-4), in which the individual changes 
[du dv dw)’^ should be replaced by values corres|K>nding to the particular problem 
at hand. Since the number of possibilities is unlimited (note that this general approach 
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may also be applied to "spat ial or three-dimensional geodesy"), in this section only 
two examples will be shown. 

4.6.1 Differential Changes in (X, 0 , h) Due to a Change 
(d Xo, d<Po, d h(l at the Datum Origin (X<„ <0o. ho) 



42 



the rotated curvilinear frame parallel to the particular local system at any point P. 

The complete operation may be expressed symbolically by the following 
commutative diagram: 

(T). i, c)o - (t?, c. C)t 



(u, V, wlo 


4.6.2 Effect of Rotations on the Curvilinear Geodetic Coordinates 
From (4.1-4) and (4.2-10) 


J— 1 


p 



(N +h)cosodX 


0 6w - 6(|) 


(N +h)cos<pcosX 

(M + h) d<p 

= R 

- 6o) 0 


(N +h)cos<psinX 

dh 

J 

6r 

6^ - 6c 0 


(N(l - e^) +h]sino 


(4.6-6) 

After the above matrix multiplications are performed and after simplification (see 
Appendix C), the following three equations are obtained: 

dX.60 = - 6co + 6€^1 - ^ tanocosX + 6({, ^1 - J^^jtaiupslnX (4.6-7a) 

d<p6H = -0€sinX ■ ' - omcoaX (4.6-7b) 

M + n M + h 

dhfin = - 6€Ne^sin<pcos(DsinX + 6({)Ne^sinocosoco8X (4.6-7c) 

The above equations in similar form are also given in (Hotine, 1969, p. 263]. b 
Appendix A the reader may find in equations (A. 1-1^ the effect of rotations as given 
in (Molodenskii et al. , 1962]. Note that equations (4.6-7) are not completely rigorous 
expressions because the assumption of small rotations is implicit in the matrix 6 U. 

See equation (4.2-10). 

The complete rigorous expressions in matrix iu>tation may be obtained if in 
place of 6R in (4.6-5), the matrix [U - I] is substituted, where R is given by (4.2-6). 
The utility of equations of the type (4.6-6) in the application of minimal constraints to 
curvilinear coordinates is discussed in Appendix D. 
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5. APPLICATIONS TO OTHER CURVILINEAR SYSTEMS 


5.1 Ellipsoidal Coordinates 

The tbeoiy developed in Chapters 3 and 4 was applied exclusively to a case of 
orthogonal curvilinear coordinates, the geodetic coordinates (X, (p, h). It will be 
shown at this time that the above theory may also be used in other orthogonal cur- 
vilinear systems of common application in geodesy and geophysics, such as ellip- 
soidal and spherical coordinates. 

Although in both instances, the family of surfaces is triply orthogonal, one 
basic difference, however, should not be overlooked. It is simply that while the 
coordinate lines and surfaces generated by geodetic coordinates are orthogonal, they 
are nevertheless not confocal. In this chapter only confocal families of surfaces will 
be treated. The following implications hold: 

confocal ==> orthogonal 

orthogonal confocal 

Appendix E reviews the properties of some common families of ellipsoids. 

For the sake of generality, it will be convenient to give first the transforma- 
tion equations between the Cartesian and general ellipsoidal coordinates (reference 
ellipsoid of three parameters) from which, as is known, two degenerated cases can 
be obtained: rotational ellipsoidal (reference ellipsoid of two parameters), also called 
spheroidal by some authors, and spherical coordinates. 

On the following pages, it will be assumed that the reference ellipsoid for the 
ellipsoidal as well as the curvilinear geodetic coordinates is geocentric. Therefore 
the notation (x, y, z) will be used for the Cartesian system. 

The transformation equations between the different curvilinear and Cartesian 
coordinates are given l)elow. It is assumed that the reader is familiar with the basic 
definitions. Otherwise (Hobson, 1931] or [Heiskanen :uid Morit/., 19(>7] can be consulted. 
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b) Local Ellipsoidal and Geodetic Frames 
Fig. 5.1 Ellipsoidal Coordinates (after [Helskanen and Moritz, 1967.]) 
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(5.1-1) 


(5.1-2) 


(5.1-3) 


The elements of the Jacobian matrix for the transformation between general 
ellipsoidal and Cartesian coordinates may be computed easily: 


— = - (u" + E^)^cosj?sinX 

bx 

(u“ + E )‘sin2fcosX 

~ = 1 cos^ccs X 

dX 

0)9 


du (u^ + E")® 

^ ~ (u^ + Ef)^cos^cosX 


(u‘ t Ei^)^sln^3inX 

if Tcosj?sInX 

dX 

d/9 


du (u‘ +Eir 

^ = 0 


u cos ^ 

~ = Bin'0 

dX 



d u 


(5.1-4) 
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and then considering equations (3.2-12) and the above results, the elements of the 
metric matrix when the reference ellipsoid has the parameters (u, E, Ei) are 
obtained as follows: 


cos'0[u^ + E'^sin^X + Ei^cos'^X]^ 

(5.1-5a) 

[iP + sin^jS (E^cos^X + Ei^sin^X)]^ 

(5.1-5b) 

bcos®?r/ 

(5.1-Gc) 

L \u“+E^ u-H-Ei*-/ J 


Assuming now E E i , the elements of the Jacobian and metric matrices for 
the transformation between ellipsoidal (rotational) and Cartesian coordinates are immed- 
iately obtained 


|~ = - ( u^ + E ^ )^cos ]8 sin X 

O A 

= ( u^ + E ^ )^cos /3 cos X 


= - (u*^ + E^)^sin^cosX 
op 

1^ = - (IP +E^)^sinj88inX 


^ = — -- — J. cos ^ cos X 

du (u“+E^)^ 


u 


1 

Bu (u^+E"")^ 


cos /3 sin X 



d Z a 

Jg - u cos S 


— = sin /9 
5u 


(5.1-6) 


and 


hi 

h2 

ha 


cos/3 (u^ + E"")^ 

(u^ + E"sin'^i3)^ 

1 

[;3— cos' j3 t 8in‘' j8]® 

u'- + E'' 


- ( 


u" + E^sin^ ff 

~3 . .,3 

u + K 


= lU 


(5.1-7a) 

(5.1-7b) 

(G.l-7c) 


Finally, for the simple case of spherical coordinates (E = 0), the following 
known relations are obtained: 
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= - rcosj^)sin\ 

O A 

1“ = - r sin(()cosX 

= cos(|!)cos X 


1“ = rcos^eosX 

0 A 

= - r sin0sinX 

1^ = cos 0 sin X 

(5.1-8) 


d z , 

jj = rcos* 

d z 

- = sin* 



h 1 = r cos i/) 


(5.1-9a) 


ha = r 


(6.1-9b) 


ha = 1 


(5.1-9C) 


5.2 Transformations Between the Normal Gravity Vector Components 

Because of the nature of the reference body, the normal gravity field is gen- 
erally represented in function of ellipsoidal coordinates. The function U(tt, P) is 
given for example by equation (2-62) in [Heiskanen and Morit^ 1067] for a particular 
rotational ellipsoid with semimajor axis a, ’ '-^ar excentricity K, gravitational 
constant kM and rotational velocity co. 

Assume now that the following transformation is desired: 

(Vx, Vy. Vz) (ycr.yf.Vs) (6.2-1) 

t) ^ C 

where 

^ ^ r- / . /SU au c)U\ 

y = gradC 

That is, given the components of the vector y (treated as a free vector) in the 
geocentric system (x, y, z), obtain the components of the normal gravity vector 
In a local ellipsoidal frame (T), O at the point P(X, )3, u). 

The above transformation will be performed through an orthogonal rotation 
matrix. Therefore the inverse transformation (which is really the one practically 
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used) is immediately available. 

To obtain the rotational matrix in the transformation (5.2-1) is not 
as simple by geometrical deductions as in the case of geodetic coordinates. 
The main reason being that it is very difficult to visualize the connection between the 
angle /3 and the local ellipsoidal frame (see Fig. 5.1). 

hi cases like this, as mentioned in Section 3.8, equations (3.8-1) and (3.8-2) 
may be applied. For this particular curvilinear system, the following relation holds: 

R, = H/j/ (5.2-2) 


where the elements in the matrices J„ and are given by equations (5.1-6) and 
(5.1-7) respectively. Thus 


y| 


~y* 


y* 

= R, 

Yr 


Yr 

7 c. 


Yt 


Yz 
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l/(u^ + E‘^)^cosj8 


- (iT® + E^)^cos/3sinX 

- (u^ + E*^ )^sin/2cosX 
u cos/3cosX/ (u^ +E^)^ 


l/iu(u'’ +E“)^ 


( u ^ + E^)KobPcob\ 

- (u^ + E"^)^lni38inX 
ucosjSsinX / (u^ + E^)^ 


(5.2-3) 


l/ui 

0 

ucos/9 

sin/3 


y* 

7y 

Tr 


(5.2-4) 


The above is also given in equation (6-11) of [Heiskanen and Moritz, 1967], although 
the final result Is obtained by a different approach. Observe that the notation of the 
equations presented here is in accordance with the general criteria of this report. 


49 



After the matr'x vi ultiplication in (5. ?-4) is performed, the rotational matrix 
of the transformation (5.2-1) is as follows: 


- sinX 


cos X 


R- = 


sinjScosX 

Ul 


— sinjSsinX 

lU 


u 


-TCOSiScOsX 




ui(u ’ + E ' ) 


ui(u" + E")^ 


cos /3 sin X 


u a 

rcosp 

uj(iP +E'-’)^ 


sin/3 


Ul 


J 

(C.2-5) 


5.3 Differential Transformations Between Cartesian and Ellipsoidal Coordinates 


The transformation between differential changes in Cartesian and ellipsoidal 
coordinates can immediately be obtained trea^iug tne basic eo-iations of Section 4. 1 in 
a more general way. 

Applying (4.1-3) to this specific case, 



dX 


dx 

II , 

d^ 

- H, 

dy 


du 

— ^ 


dz 


but recalling (3.8-1) and (3.8-2) 

U, = (5.3-") 


one obtains 
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(5.3-3) 


dX 


dx 

d)3 

= (h;')=^j/ 

dy 

du 

_ 


dz 


dX 


dx 


diS 

= 

dy 

and ^ 

du 


dz J 



dX 


dx 

di9 

= J,' 

dy 

du 

j 


dz 


(5.3-4) 


Observe that equation (5.3-4) was expected, if the total differential approach is recalled. 
From (5.3-3) the analytical form of the inverse of the Jacobian may be computed, 

giving 


- sin X/cos^ 


ui^(u'^ + E"^)' 


- sin /3 cos X 
ucosjScosX 


ui^cosX/cos/3 


- sin)3sinX 


ucos/SsinX 


0 

ucosff 

• > D 

(u' + E r 
slni3(u" + I? )^ 
(5.3-5) 


5.4 General Commutative Diagram for the Transformation of Free Vector Components 

Figure 5.2 represents all possible transformations or mappings between free 
ve * components in spatial rectangular coordinate systems. Emphasis has been 
placed on the coordinate ^sterns discussed in this report (mainly geodetic and rota- 
tional ellipsoidal) but the same logic n'.ay be applied to any other orthogonal curvilinear 
coordinate system. 

For an easy recall, some of the basic matrices represented in the diagram are 
given according to equation number in the following table; 
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Table 5.1 Equation Numbers of Commutative Diagram Matrices 


Matrix Type 

Geodetic System 

Rot. Ellipsoidal System 

i^mbol 

Equation # 

symbol 

Equation # 

Jacobian 

j 

(3.2-8) 

J, 

(5.1-6) 

Jacobian Inverse 

j" 

(3.7-2) 

J.' 

(5.3-5) 

M^ric (Diagonal) 

H 

(3.2-12) 

H. 

(5.1-7) 

Rotation (Orthogonal) 

R 

(3.3-2) 

R, 

(5.2-5) 


As a simple application from the diagram, let's assume that the com- 

ixjnents (y , y?, y;y)are desired. The diagram gives immediately 
T? ^ C 


3F / dq 


dF / bX 

dF / di 


bF / bP 

dF / dC 


bF / bu 


where F is any scalar function. 

Assumii^ F to be the spheropotential function U, one can obtain the 
components of the normal gravity y along the local ellipsoidal system at the point 
P(X, |8, u) 


y.v 


eu / dX 

V 





dU / d/3 

M. 


dU / du 
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Thus the \ 'ues of (y ^ , y j-, y^) can be computed readily from the metric matrix of 

1? ^ C 

the transformation between ellipsoidal and Cartesian coordinates and the partials of 
the known function U (u , /3) respect to the ellipsoidal coordinates. Clearly, for a 
rotational ellipsoid BU / dX = 0. 

Once y (y~, y?, yj) is computed, it is immediate to obtain (y*, y,, yi) as 
s C 

explained in Section 3.2, through the transformation 




It' 


(yK * » yt ) 


which is also implied in Fig. 5.2. 

Observe that the diagram of Fig. 5.2 is commutative. This means that there 
are several ways to obtain the transformation between two sets of free vector compon- 
ents. The selection of the approach will depend on the type of data readily available. 
For example, the differential transformation between ellipsoidal and geodetic coordi- 
nates may be performed according to the following possibilities. 


(dX, d/3, du) 
(dX, dj3, du) 
(dX, d/3, du) 
(dX , d^, du) 


H 'Ult ^ H, 


- .a, ■ -fl. - > 




j'j. 


(dX , dv), dh) 
(dX, d<p, dh) 
(dX, d<p, dh) 
(dX, d<p, dh) 


where the matrix of the transformation is given by: 


J^J 


((u ' + E^)^sin^sin<o 
+ t5rcos/3cos<p] / (M -t h) 


- (iP + E' ) ^sin3cos<p 
+ uT cosj3sin(P 


0 


t-u(u"+E"’) ^cos)3sln<p 


+ sinj3cos(Pl / (M + h) 

■ - i 

u (u" + E") cospslno 
+ sin^sin<p 

(5.4-3) 

Notice that there is no change in the transformation of dX between ellipsoidal and 
geodetic coordinates, as expected. 
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5.5 The Rotation Matrix ^ Between the Geodetic and ElHpaoirtfll Local Systems 


At this point it will be interesting to stuc(y the nsformation between the local 
geodetic and ellipsoidal ^sterns. That is, we are interested in obtaining the rotation 
matrix ‘’R of the transformation 


(n. (ti.C.C) 


From the diagram o'J Fig. 5.2, 



Thus, after matrix multiplication. 




RR/ 



0 
m 
- n 


0 

n 

m 


where 


m =i[u(u"+E^) ^ cos |3 COSO + siniSsittip) 


\ i 

n = — (-u(u +E") cos^sino + siniScosoI 


(5.5-1) 


(5.5-2a) 


(5.5-2b) 


(5.5-2C) 


The matrix It will transform components of the normal gravity vector y from the 
local ellipsoidal to the local geodetic system, or viceversa. The orthogonality of 
may be proved easily. 
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Notice that represents a counterclockwise rotation a along the 17 axis (Fig. 5.1b) 
where 


cos a = m 


and 


sin a = n 


This can also be shown by an independent approach considerii^ the orthonormal bases 
(© 1 * ea. 63 ; and (Ci, 63 , 03 ) along the respective local frames. 

Clearly , 


cos Of 


-> -♦ 


Ca . §a 



-» -» 


C3 • Oa 



The components of the vectors of the base e 1 (i = 1, 2, 3) may be found in 
equations (3.2-13). They are also given by the row elements of the matrix It in 
equation (3.3-2); the components of the vectors in the base e j (i = 1, 2, 3) are the 
elements of the rows of ll^ given by (5.2-5). 

Considering that these bases are orthonormal, that is, all the vectors are 
unit vectors, it follows that 


COSO! 


es 


= e: 


- [u(u 

UJ ^ ' 


+ E ' ) ^ cos /3 cos o + sin /3 sinol 


Therefore , after algebraic manipulation and simplification, 

sin‘'a = 1 - cos' Of = ^ [ (u~ + E")^sin/3coso - ucosjSsincp] ’ 

u= +E^sin^i8 


Thus, finally 



5.6 Transformation between Local Geodetic and Local Astronomic Systems 


The same criterion of the above section may also be ^plied to obtaining 
the transformation of vector components between the local geod^ic (t), 4 > C) 
and local astronomic ^ Z*) systems. 

Clearly in this case, the mapping 


(n. C. O > (n*. i*, C*) 


will be obtained as follows: 


where 


and 




n* 


n 

Z* 

= R*R^ 


cv 


.'J 


(5.6-1) 


R* = Ri(90 -<p*)R3(90 +X*) 

<p * = reduced astronomic latitude (5.6-2) 

X * = reduced astronomic longitude 


cos(X* - X) 

- sin<p%in( X* - X) 

cos<p*sin(X* - X) 


sin^sln(X’" - X) 

C08<pC0S<p* 

+ sin^8in<p*cos(X* - X) 
cos<psin<p* 

- sin<pcos(p*cos(X* - X) 


- cospsin(X* - X) 
sin<pco8<p 

- COSO) 8ln(p*oo8(X* - X) 

sin(0 sin<p * 

+ cos 0COB<P * C 08 ( X * - X ) 
(5.6-3) 


Assuming now small differences between the geodetic and astronomic coordinates: 


(p* - (p + 6<p (5.6-3a) 

\* = X + 6X (5.6-3b) 


and with the simplifications 
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8in(X * - X) fta 6 X 
COS( X* - X ) sa 1 


(5.G-4) 


8in(<p*-<p) at 6<p 

COS(^* - <£>) Ri 1 


COSO" 


CO8<0 


it follows 


1 siac>6X - coa<p6X 

RR*^ fis - sino6X 1 - 6<p = Re* (5.6-5) 

cos<p6X 5<p 1 

This is really the case when one wants to transform between the local geodesic 
(V t ^ . C) and local astronomic (tj*, 4 *» C ’*' ) systems at some iiarticuiar station 
P where the deflections of the vertical T) ' and 4 ' are known. 

Then, 

4' = 0* -<p ~ 6(p (5.6-6a 

r\' ~ (X*-X)cos<D = 6Xcos<p — ■ sin<p6X ^ Tj'tano (f>,6-6li; 


Therefore, substituting in (5.6-5) 


r\ ' tan«5 


RR*' 


T)*tan<p 


-e* = Rft' 


(5.6-7) 


Finally, if differential changes in (r), 4 t C) due to small rotations of the local geodetic 
system by amounts ( 6X, bo) at the point ( X, o) are sought. 


d4 = [Rft* -n 4 


where Re* is given by (5.6-5), thus 


dt) 0 

d4 - - slnoSX 

d 4 cos o 6 X 


sin<p6X -cos<p6xl T t) 


(5.6-8) 


4 1= 6H* 4 


(5.6-9) 



6. SUMMARY 


The present work uses a general matrix approach in reviewing some basic 
differential transformations between Cartesian and curvilinear coordinate systemd. 

The md^hods discussed here are applicable to any type of orthogonal curvi- 
linear coordinates. Nevertheless in this report only geodetic and ellipsoidal (rota- 
tional) coordinates are examined. 

As an application of the theory, differential changes in geodetic coordinates 
due to shift, rotation and scale of the geodetic system are found. 'I'he same results 
may be obtained employing other methods, such as the total differential approach or 
tensor calculus. 

This study also tries to clarify the confusion in recent geodetic literature re- 
garding the so-called "Molodenskli model," which is used in the least squares solution 
of the seven transformation parameters between world and geodetic (datum) systems. 
Careful consideration of the differential equations given in (Molodenskli et al. , 1962] 
shows that the model attributed to them is not impliv . ’r. i .eir work. 

After defining three basic transformation matrices (rotation R, metric H and 
Jacobian J), mappings between differential changes in the Cartesian jind orthogonal cur- 
vilinear coordinates are established. This is illustrated by a general commutative dia- 
gram (Fig. 5.2). As an example, the following differential transformation is presented 

(dX, d/3, du) < ^"(d X , do, dh) 

(To the best knowled^^ .'•f author, the transformation matrix of this mapping is 
given here for the fi. st t,‘«r j) Thus, it is possible to obtain differential changes In 
the geodetic coordioaced ct a point P(X, o, h) P(X, /3, u), when the ellipsoidal 
coordinates change differentially or viceversa. 

As an immediate result of the above-mentioned diagram, it is possible to obtain 
the components of the normal gravity vector in the local geodetic system; thus the com- 
ponent of attraction along the geodetic normal can be found. 

Finally, the rotation matrix between the local geodetic and local ellipsoidal 
frames is examined. 
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APPENDIX A (Referenced in Section 4.4.2) 


A. I Matrix Form of the Equations Given in fMolodenskii et al. . 19621 

Those acquainted with the English translation of the work [Molodenskii 
et al. , 1960 ] know that the approach followed there uses the strict procedure of 
differentiation of curvilinear coordinates on practically the same line as was intro- 
duced in the classical work by Lam^ [1837]. Therefore, no mention of local or mov- 
ing frames is evident in the Russian translation. Thus, in order to change over to 
the matrix notation of this report, a correspondence between frames and their rota- 
tion matrix must be established. This is shown in Fig. A.l. 


z 



Fig. A. 1 Local (loodilic Frame in Moloder ’.ii 



Thus the rotation matrix OR may be written: 

cos<Oc*osX cos(pslnX sincp 

IR = Ri(-<p) Ra(X) = - sinX cosX 0 (A.l-l) 

- sin<pcosX - sinipsinX cos<p 

In the equations given in [Molodenskii et al. , 1962] the differential changes (da, dbi 
in the semi-major and semi-minor axes of the reference ellipsoid are introduced. 
However, in Section 4.4.2 of this report the flattening was used instead of b, thus 
the following substitutions must be taken into consideration 

f = => b = a-r*f (A. 1-2) 

a 

and differentiating above 

db=da-adf-fda db=da(l f)-adf (A.l-Jl) 


Therefore 


and finally 



(A. 1—1) 


The exact correspondence between equation (1.3.2) in [Molodenskii et aL, 1962] and the 
matrix notation used in this report is as follows (see also Eq. 4.4.5): 







3x 

dx 

aa 

d f 



da 

d f 

dz 

dz 

da 

d f 



(A. 1-5) 
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where now from (4.4-3), making use of (4.4-4), (A. 1-2) and 



1 • e = T 
a 


' possible to write 



1) 

* M8in^0codOsinX 

(A.l-6a) 

(1 

§ Msin^cpcososinX 
b 

(A.l-6b) 

11 

b a 

_ (M - aN - Mcos o) sino 
a 

(A.1-6C) 

hi accordance with the notation in [Molodenskii et al. , 

1962] , the following 

equalities are established. 

P = aW 

(A. 1-7) 


N = 5^' 
p 

(A. 1-Ha) 


II 

(A.l-Hb) 


Premultiplying both sides * equation (A.. 1-5) by the rotation m;itrlx lit and 
recalling (4.4-8), equation (1.3.3) in [Molodenskii et al., 1962] follows immediately. 



dx 


dx 


IK 

dy 

= K 

dy 

+ R 


d z 

p. — 


dz 

1 

L 


0 

- fz 


0 

- c. 


- fy 
f, 
0 


X 

y 

■ 1 . 


dh 1 
(N + h)eoS(pdX 
(M h)do 


+ IK 



dx 

da 

d f 


dy 

da 

d f 

d/. 

dz 


d f 


da 

b / ^ db \ 
a \ a ^ ' 


(A. 1-9) 
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It can be r 1 that Equation (A. 1-9) corresponds to (1.3.3) given in (Molodenskii et al,1962] 
Fc c mple, it is immediate to show that the terms corresponding to the rotations 
in (. 9) are equal to the ones presented by Molodenskii et al. , (1962). 


(1) ' 


0 Cr - €y 


X 

(2) 

= IR 

- 0 Cx 


y 

(3) 


€y — Cx 0 


z 


(A. 1-10) 


costpcos X 

cos to sinX 

sintp 


0 


- €y 


X 

- SinX 

cos X 

0 


- Cx 

0 



y 

- sintpcosX 

- sintOsinX 

cos to 



- 

0 


z 


- costOsinXCj + sintOCy 

cos0cosX€i - sinto^K 

- costpcosXCy +co8<pslnXCx 

X 

- cos Xe. 

- sinX^x 

sinXCy +cosXCx 

y 

sintOsinXCi +cos0€y 

_ 

- sin to cos X^x - cos to fx 

sin0cosXfy - sinosinXe. 

z 

— 


X 

y 


(N + h) COS0COS > 
(N +h)cos 0 sinX 


(N +h)cos0cosX 
(N +h)cos0sinX 

z 


[N(l - e^ ) + h] sin0 


^ + h^ sin0 


thus 

a) 


= -(N + h) cos^<psinXcos XCt + (N + h) sin<D cos tocos Xc + (N + h) cos^tOsinXcosXCi 
-(N + h) slop cos tosinXf, ^ + h jsintocostocos XCy 

+ N + l^s into cos (psinXfx 

= AsintOcostO(»'yCosX - f, sinX) 


where 

A 


= (N +h) N + hj = n/i - ^ N ^ ^ 


a‘ - b^ 


therefore 


( 1 ) 


a - b“ 


F’ 


sintOcostO (CyCosX- CySinX) 


(A. 1-11 a) 


06 



Now 


(2) “ - (N +h)cos<pcos“Xc, - (N + h) cososin^Xe, ^ sin<psipX< 


K N +h^ si 


sinocos X€» 


= - (N + h)cos<p€j+( ^ N + hj sin<o(€^xCOsX + CySinX) 


and 


(2) = - (N +h)ooso€i + z(€yCosX + fySinX) (A.l-llb) 

Finally, 

(3) = (N + h) sinoc ^^^sinXcos Xf, + (N +h)cos^OcosX€y 

- (N + h) sinOf;os<psinXcosX Cj - (N + h) cos'^OSinXc* 

+ [N(l - e^) + h] sin^ ocos XCy - [N(l - e^) + h] sin“ (OsinXc, 

- B(€j.cosX- e^sinX) 

where 

B - (N + h) cos“0 + [N(l - e’) + h] sin“ o = N+h-Ne" sin^ O = N + h = p + b 
Thus, 


(3) = (p +h)(€yCOsX - c^sinX) 


(A.l-llc) 


To conclude, one may write the effect of differential rotations c,, €y , on the 
geodetic ooordinates as given by Molodenskii et al-, [19G2J 


(A.l-12a) 
(A.l-I2b) 
(A.1-I2C) 

These equations are equivalent to (i.G-7) with 

e, 6 <■ 6C (r 

and taking into account (A. 1-7) anc' (A.l-S) 


(N +h) cosod X = - (N + h)cosof: + z(e,cosX + e^sinX) 


(M+h)do = (p + h) (e.cos X - e,sin X) 
a - iy 


dh 


sinocosolfycos X - e^sinX) 
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APPENDIX B (Referenced in Section 4.5) 


B.l Differential Changes in (X« <p, h) Due to a Change 6L in Scale 


From (4.5-1) 


(N + h) cos<pdX 


- sinX cosX 0 


(N + h)cosocosX6L 

(M + h) do 


- sinocosX - sinosinX coso 


(N +h)cososinX6L 

dh 

1 6l 

cosocosX cosOsinX sino 

L j 


(N(l - e ) + h] sinofi L 


After matrix multiplication 

(N -- h) cos<pd Xsl = - (N +h)cos<osinXcosX6L + (N +h)co8<psinXcosX6L = 0 
Thus, 


(M + h) d<D6L 


Therefore, 


dX6t = 0 (1.5-2U) 

= - (N + h) sin<pcos<ocos‘ X6L - (N +h)siaoco8<psin^X6L 

+ [N(l - e") + h] 8in<pcos(p6L 

= - (N + h) 8in<pco806 L + (N +h)sin<pcos<p6L - Ne^8in<pcos<p6L 
= - Ne‘^sin<ocos 06 L 


d06L = - 


Ne^sin0co8<o 
M + h 


6L 


(4.5-2b) 


Finally, 

dh = (N + hlcod^cDcos"^ X6 L + (N + h) co8“<psin“X6 L + [N(l - e^) + . in‘*<p6L 
= (N +h)cos“<p6L + (N +h) 8 in- 06 L - Ne'^sin~<£)6L 
= (N +h)6L - Ne’sin'^oSL 
= (N(l - e^ si.i^ O) +h]6L 


but making U8e of equations (4.4-4) it follows that 


dhst = (aW +h)6L 


(4.5-2C) 
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B.2 Differential Changes in (X, o, h) Due to a Chantie 6a in the Semimalor Axis 
of the Ellipsoid 


From (4.5-3) 


(N +h)cos0dX 


- sinX 

cosX 

0 


COS0COSX - 

w 

(M + h) d0 

= - 

- sinocosX 

- sin0sinX 

COS0 


co8<p8lnX ^ 

— w 

db 

5 a 

COS0COSX 

cos08inX 

sin 0 


(I - e^siixp 

L w 


, ,, co8<oainXcosX . cosipsinXcosX . . 

(N +h)cos( 0 dX 5 , = 6a 6a = 0 

w w 


Thus , 


(M +h)d(0 6. 


Therefore* 

Finally, 

dh5, = 


dXft, = 0 


(4.5-6a) 


sin<pcos<pcos^X , ^sinocoscpsiri^X , (1 - e^)sin<pco8(/:> . 

w wf w 

sin(pcos<p . (1 - e^) sin(pcos 0 . 

w 

e^sin<pcos(P , Ne^ sin<pcos 0 . 

5a = 6a 

W a 


, Ne“sin0cos0* 

d0 6* = - ,v ;- 6a 

(M +h)a 


(4.5-Gb) 


cos^pcos'^X . cos“0sin®X- (l-e'lsin'^ . 

— 6a rr: 6a m 6a 


W W 

cos^o , (1 - e')sin‘^(0 


W 


w 

1 - e^sin^0 
W 


W 


6a 


6a 


Thus, 


dhft, ’-VV6a 


"-Gc) 
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APPENDIX C (Referenced in Section 4 . 6 ) 


C.l Differeatial Changes in (X. <p. h) Due to Small Rotations 6 f. 6 iL 6 u} 


From ( 4 . 6 - 5 ): 


(N + h) cos(( 3 dX 


0 6(0 - 6(/j 


- 

(N + h) cos (p cos X 

(M + h)d0 

= U 

- 6(0 0 6 < 


(N ^h)cos(PsinX 

dh 

6 R 

6 i() - 6 € 0 J 


[N(l - e ) + h] sin<p 


- cosX 6(0 

sin<osinX6(o 
+ cos( 06 V> 

- cos(£)smX6co 
+ sin^6((i 


- sinX6co 

- sii ocos X6(0 

- cos(p6r 

coscocos X6(o 

- sinc)6€ 


sinX6^ + cos X6c 

sin0cosX6(J) 

- sin<DsinX6 € 

- cosOcos \ 6 il) 
cos<psinX6 C 


(N + h)cosOcosX 
(N + b) eos<psinX 
[N(l - e ) + hj sinO 


(N +h)cos<odX6R - - (N + h)cos<pcos' X6(o - (N + h) cos< 5 sin“X 6 <o 

+ sinXsin<o 60 (N(l - c") + h] + siiuOcos X6c [N(l - e~) + h) 


= - (N +h)cos(06(0 I [(N + h) - Ne ]sino(sinX 60 ^cosX6e) 


Then 


.X <■ (N + h) sm<p - N e sin<0, . . * , ^ x * . 

dX6R = -6(0 + (sinX6^t + cosX6€) 

(N + h)cos(£> 


= - 6 (0 + tanc 3 ( 1 


Ne^ 


N +h / 


\ (sin X6(() + cos X6 f ) 


Thus, finally: 


/ N e \ / N e" \ 

dXjR = - 6(0 + 6e (1 - |tan<ocosX + 6||; (1 - 1 tan0sinX ( 4 .G- 7 a) 


(M +h)do6R (N +h)cos<ocosXsln(OsinX6(0 + (N + h)cos' <pcosX6(i) 

- (N + h) cos(P8inXsin<PcosX 6(0 (N ' h) cos“ 08 tnX 6 f 
+ (N(l - c“) + h] 8in <^cosX6(() - (N(l - o") + h] ain“ 08 inX 6 f 
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(M +h)d<p6R = (N +h)cos^<pcosX6^() + (N + h) sinP<pcosX60 - Ne^ sinP<0cosX6^ 
- (N + h)cos^V>sinX6 € - (N + h) sia?;c>sinX6 f + Ne^sin'<p8inX6€ 


= 6^^)cosX((N +h) - Ne^ sin'(p] - 6csinX[(N +h) - Ne^ sin^<p] 


But, 

(N +h) - Ne"^siiT0 = N(1 - e'^sin^ip) + h = NW® +h = +h = aW +h 

w 

Therefore, 


(M +h)d0 = 6^cosX(aW +h) - 6€sinX(aW +h) 


and finally. 


d<p6R 


- - 6 f sin X 


aW + h 
M +h 


+ 6 ill cos X 


aW +h 
M +h 


(4.6-7b) 


dhfiR = - (N +h)cos^</jcosXsinX6w + (N +h)cos<pcosXsin<p6^() 

+ (N +h)cos^<psinXcosX6t»j- (N +h)cos(0sinXsin<p6e 
- [N(l - e^) + h) sin0cos<pcos X60 + [N(l - e^) + h] sin<pco8<psinX6€ 


= - 6l(i sin<pcos<pcosX[(N +h) - Ne^ - (N +h)] 
+ 6c sin0cos<osinX[(N +h) - Ne“ - (N +h)] 


Thus: 


dhjR = - 6cNe^sin<0cos<psinX + 60Ne^sin0co80cosX (4.6-7c) 
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APPENDIX D (Referenced in Section 4.6.2) 


D.l Innpr ^Minimal) Constraints in Curvilinear Coordinates 

D.1.1 Introduction 

Papers on minimal constraints and their application to ^odesy are abundant 
in the literature. The basic principles introduced here follow [Pope, 1971] where the 
interested reader can consult the fundamental references on this topic. 

It is well known that in most geodetic problems the set of normal equations 

NX +U = 0 (D.1-1) 

is a singular system when the original observation equations F(X, L) = 0 do not con- 
tain sc me peculiar constraints. 

In the specific case of a spatial network the following relations hold 

nNn => Rank(N) = r (D.1-2) 

where r < n and e = n - r - 7 

The value e is generally called the rank deficiency of N or the degrees of 
freedom of the network (not to be confused with the concept of degrees of freedom 
a least squares adjustment [sec Uotila, 1967]). 

As a consequence of (D.1-2) 

{ N 1 = 0 ==*> N is singular 

One way to solve equation (D.1-1) in this case is by bordering the normal matrix 
N and solving the system : 


' N 

E ' 


X ■ 


' - u ■ 

. E' 

0 


_- K _ 


0 

L 


where n^e is basis for the solution space (null space of N) of the homogeneous 
equation. Therefore 

NK--0 (D.1-4) 

Obviously property (D.1-4) also implies 

AE -- 0 (D.1-5) 
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The complete solution of the normal equation (D.1-1) in the case of singularity is: 

X = X, + E^ (D.1-6) 

where 


Xo 

and j3, 

NX = 0. 


any particular solution of NXo = - U 

complete solution of the corresponding homogeneous equation 


Minimal constraints are the smallest number of constraints e that produce a 
nonsingular matrix M, and minimizes X ^ X, 


M = 



(D.1-7) 


Of all the minimal constraints possk'c, some have simple geometric interpretations; 
this subset of minimal constraints is called "inner constraints " [Blaha, 1971]. 

As an illustration, assume that in E^ only angles are measured In order to 
establish a network of points. Clearly, the degrees of freedom of this network will 
be seven, if one considers that translations, rotations and scale variations will change 
the coordinates of the points, although without affecting the values of the measured 
angles. In other words, one may sa. that coordinates are not estimable quantities. 
Thus, in this example 

e = T + P. ♦•8 = 7 


where 

T = number of constraints required for origin = :i 

It number of constraints required for orientation :$ 
S number of constraints required for scale - 1 

Therefore, 

~ El Ej) nx7 

nx3 nx.'^nxl 


(13.1-8) 


D.1.2 Inner Constraints in llect angular Coordinates 


Still following [Pope, 1971] the set of inner constraints when rectangular coordi- 
nates are used, may be obtained through the differential changes in the ('artesian 


7;t 



coordinates due to translations (shifts) 


rotations and scale, that is. 


dx 

dy 

dz 


Ax 

Ay 

Az 


dx 

dy 

dz 


0 

6co 

- 60 

- 6co 

0 


6^ 

- 

0 


dx 


X 

dy 

= 6L 

y 

dz 

s 

z 


(D. l-9a) 


(D. l-9b) 


(D.1-9C) 


Thus, one can write; 


X 




Ax' 


X 


Xo 


Ax 

y 


y 0 


Ay 

==> 

y 

= 

Yo 

+ I 

3 x3 

Ay 

z 

1 

- 

1 

Az 


J 

t 

Zq 

1 

Az 


(D.1-10) 


which is in the form of (D.1-6) and gives. 


1 

0 

0 


Ax 

0 

1 

0 

; i3 = 

Ay 

0 

0 

1 

1 

Az 


In the same way, the inner orientation constraints can be found. 


X 


Xo 


0 - z y 


■ 6 (~ 

y 


yo 


z 0 - X 

1 

60 

z 

u 

1 

_Zo _ 

1 

- y X 0 

1 

6co 


(D.1-11) 


(D.1-12) 
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Therefore , 


0 

“ z 

y 


~6c 

z 

0 

- X 

; = 

6ili 

-y 

L. 

X 

0 

_ 

1 

6(A) 


For the inner scale constraint one has 


X 


Xo 


X 

y 


y<> 


y 

z 

i- J 

1 


1 

z 


Thus, 



and /3 = 6L 


(D.1-13) 


(U.1-14) 


(D.1-15) 


Finally the matrix (D.1-8) is given by 



1 

0 

0 

1 " 

- Z 

y 

i X 


El = 

0 

1 

0 

1 Z 

0 

- X 

1 

1 V 

(D.M6) 

3x7 




1 






0 

0 

1 

1 

X 

0 

1 z 

1 


D.1.3 Inner Constraints in Curvilinear coordinates (Si)herical Case) 


Curvilinear coordinates are always referred to some basic surface which 
introduces restrictions in the number of degrees of freedom needed for solving the 
network singularity. 

For example, in the case of a flat surface (plane), the degrees of freedom of 
an angular network are only four, 

Cp = T + it + S - 2+1+1 - 4 
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This is also true for spherical networks* 


63 =4 

Assuming that the scale of the spherical network is fixed through the radius r 
of the sphere* only three constraints are needed to resolve the singularity. 

Clearly the degrees (f freedom in this instance are three rotations along the x, y, z 
axis. Thus, applying (D.l-9b) to the case of spherical coordinates, one can write 
the following from (4.6-G) with 

h = 0, M = N = r, u«x, v=y and w = z ; 


or 


rcos<DdX 


0 6co - 6)^ 


X 

rd X 

= H 

- 6u) 0 6 c 

1 


y 

dh 


1 

o» 

1 

0 * 

0 

1 _ 

i 

z 


rcospdX 


0 - z y 


-6C- 

rd</3 

= R 

z 0 - X 


6 $ 

dh 

L. J 


- y X 0 


6co 


(D.1-17) 


(D.1-18) 


Therefore , 



0 

- Z 

y 


~6c' 

Er - R 

z 

0 

- X 

and ^ =• 



- y 

X 

0 


6co 


Knowing that 


'x' 


rcosOcos X 


= 

rcos<psi-. X 

z 1 

j 

r sino 


(D.1-19) 


(D.1-20) 
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Thus finally, 


rcos Xsinta 
- r sin X 


r sinXsino 
r cos X 
0 


E R 



- r cos o 
0 
0 


(I). 1-21) 


Assuming that the parameters in tbs normal equations are given by the vector matrix 

X = [dXi d!pi dXg d<Ps ... dXi d<pi... dX„ d0nT (D.1-22) 

The set of submatrices Erj required in order to avoid singularity are given in the 
following form: 


E 


? x3 


COS Xtanp 
- sinX 


sin Xtan(P 
COS X 



(D.1-23) 


The use of the submatrices Erj in the bordering of the normal matrix N for the 
solution of the singular system can be interpreted geometrically as in the rectangular 
case. It will give the ’’best" orientation to the spherical triang<jlation with points (X, <o)j . 

When only local networks on a sphere are involved, it will be more appropriate 
to rotate about a geocentric Cartesian system parallel to the local frame (t), 4 » Oo at 
the center of the network. In this case the following transformation applies. 


6e 


1 

O 

o 

i 


= Ro' 

Ho 

6co 


o 

3 

o 

1 


(D.1-24) 


and after substitution of the above in (D.1-18) one has the matrix equation. 


r cosod X 


0 

- z 

y 


6«o 

r d<p 

= U 

z 

0 

- X 

Ho’ 

Ho 

(ih J 


- y 

X 

0 


6cq, 
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Ero = R 

0 

z 

- X 
0 

1 

>. X 
1 

c 

)) 

Er R/ ; = 

o 

600 


- y 

X 

0 


1 

1 

o 

3 

*o 

j 


(D. 1-26) 


That is, 


E 


Ro 


r sinosin(X - Xo) 

- r(sinosin<PoCOs(X - Xo) 

r(sin<acosOoCos(X - Xo) 


+ COS(PCOS0o) 

- sin(0oUOS(P) 

rcos(X - X o) 

r sino„sin( X - Xj) 

- r co8<poSin( X - Xo' 

0 

0 

0 


(U.1-27) 


and finally, when the parameters are given in the form of (D. 1-22) the matrix 
„ E, (e = 3) is composed of the following submatrices: 


sin( X - Xo)tan<o 


O I 


3x3 


cos( X - Xo) 


- tan<osin<poCos(X - Xo) tanocos(PoCos(X - Xo) 
+ COS0O - sincpo 


sin(DoSin( X - Xo) 


- cos0o8in(X - Xo) 


(D.I-28) 
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APPENDIX E (Iteferenced in Section 5. 1) 


E.l Families of Rotational EUinsnirig 


E.1.1 Confocal, Similar and Quasi-Parallel Ellipsoids 


The following relations e immediately obtainable: 


E = ^/ a" - b® 



a 


E 

e - 

a 








_ ada - bdb 

(E.1-1) 

- b" 

_ h / da db \ 

(E.1-2) 

a a b / 

b(b da- adb) 

(E.1-3) 

a^'a" - h= 



From the above basic relations, it i - Possible to define the following types of ellipsoids: 

Confocal Ellipsoids . A family of ellipsoids is called confocal if 

E= constant ==> dE = 0 (E.1-4) 

From (E.1-1) the condition for copfocality is found immediately, 

, , JL ^ da db 

ada = bdb => -r* = — (E.1-5) 

b a 


or 

db = ^da (E.1-6) 

Substituting property (E.1-5) ia (E.1-2) and (E.1-3) one has. 


dfc 




(F.1-7) 


and 
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de 


da 


(E.1-8) 


E 
a' 

which is also obvious from the differentiation of e = — for E = con^ant. 

a 

Similar Ellipsoids . The .lame "similar" is applied to a family of ellipsoids 

when 

= 0 

From (E.1-2) or (E.1-3) the similarity condition follows immediately: 

1 . j da db 

bda = adb => — = — (E.1-11) 

a b 

and therefore, 

db = -dd (E.1-12) 

a 

Substituting the above equation in (E.1-2) and (E.1-3) one has 

dEa - eda =v/Tf-f da (E.1-13) 

Obviously, the same result is obtained by differentiation E = ea with e = const. 

Quasi- Parallel Elh^^oids . A family of ellipsoids is called "quasi-parallel" 

(the author v/as unable to find anywhere in mathematical literature a name for this 
family) if the following property liolds: 

da = db (E.1-I4) 

This implies: 

fpUp = a-b = constant (E.1-15) 

That is, tor any family of quasi -parallel ellipsoids, the product of its flattening by 


‘ L 
I 


d ; 

constant => ^ , 

de 


(E.1-9) 

(E.1-10) 
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its semimajor axis is constant. The value of the constant is the difference between 
the semimajor and semiminor axes of any ellipsoid in rhe family. 

From (E.1-1), (E.1-2) and (E.1-3) it is possible to obtain 

dEp = - da = — da (E.1-16) 

e c. 

j. b “ a . f . df d a . inv 

dfp = —2 — da = - - da => 7 -= (E.1-17) 

a a fa 

dep = - ^ ~ L da (E. 1.-18) 

a V a + b 

The following inequalities hold: 

I dbg I < ! dbp I <! dbc | (E.1-19) 

I dE p ! < I dEg 1 

; dfp ! < I dfc 1 

i dCp 1 < ! dec 1 

E.1.2 The Variation dh of the ueodetic Height 

After matrix multiplication equation (4.5-4) gives: 

dh = -Wda + -ll~ sin=’(pdf (E.1-20) 

w 

Therefore, the variation dh of the geodetic height, according to t>e different 
cases mentioned in the previous section, may be obtained, 
a) Confocal case : d E =0 

Substituting tne value of df given by equation (E.1-7) in (b.l-^iO), after 
simplication one has: 

(inr=-‘^? (E.l-31) 


Hi 



b) Similar case : df = 0 
Obviously, from (E.1-2C^: 

dhs = - Wda 

c) Quasi-parallel case : da = db 

This case, although more involved, is also easy to obtain: 

dhp = - Wda + ^'^^ -sin^( 0 da 
or 

a-(a-b)sin^O 

dhp = - ; ■: - — ■ ■■ - da 
ya^-(a" - b“)sin°o 


(E.1-22) 


(E.1-23) 


As a consequence of (E.1-23) one concludes that the variation of h in the case 
da = db is not constant. This is the primary reason for the name "quasi- parallel" 
for this family of ellipsoids. 

From (E.1-23 Hhe maximum value of dhp is obtained at 


where 


0 


arc sin 


V 


a 

a +b 


d hp M 


2N/aib 
a + b 


da 


(E.l - 24) 


(E.1-25) 
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